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One-Sentence Summary: Low-field nuclear magnetic singlet states allow the 
storage of nuclear spin order far much longer than the conventional 
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Abstract 

Nuclear spin singlet states may be used to store nuclear spin order in a room 
temperature liquid for a time much longer than the spin-lattice relaxation 
time constant T x . The nuclear spin singlets are unaffected by intramolecu- 
lar dipole-dipole relaxation, which is generally the predominant relaxation 
mechanism* We demonstrate storage of nuclear spin order for more than ten 
times longer than the measured value of T*. This pheuonemon may facilitate 
the development of nuclear spin hyperpolarization methods, and may allow 
the study of motional processes which occur too slowly for existing NMR 
techniques. 



One of the principal factors contributing to the wfde-r anging utility of nuclear 
magnetic resonance (NMR) is the relatively long lifetime of nuclear spin 
order. This long lifetime allows NMR spectroscopies to follow processes 
such as diffusion, flow and alow molecular motion (1,8), and to contemplate 
using NMR for quantum computation The relaxation of the nuclear 

spins back to thermal equilibrium, and hence the lifetime of the nucl^ax spin 
memory, is characterized fay a time constant Ti known as the Longit udinal 
relaxation time constant, or as the spin-lattice relaxation time constant. In 
most cases, this parameter is of the order of seconds, although it may be much 
longer in special cases. Major changes in Ti can usually only be achieved by 
changing the isotopic composition of the sample, or by changing its physical 
State. The sophisticated manipulation of nuclear spin interactions by radio- 
frequency pulse sequences, so prevalent in modern NMR (1 ), have always 
been restricted to a lame-frame set by the intrinsic T± liniit. Experiments 
which seek to enhance NMR signals by hypcrpolarization phenomena (5-9) 
have great potential for extending the range of NMR even further, but are 
always hampered by the need to use the hyperpolarized substance within a 
time of the order of T, . 

In this report, we demonstrate, for the first time, the storage of nuclear 
spin order for a time much longeT than 2\. So Far, we have demonstrated 
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this phenomenon on molecules containing pairs of protons, in which the main 
mechanism for nuclear spin relaxation is the modulation of the intramolecular 
dipole-dipole coupling by rcoleculax motion- The extension of nuclear spin 
memory beyond the Ti limit is achieved by creating nuclear spin singlet 
states in low magnetic field (20 mT or less). As shown below, these law- 
field singlet states (LFS) are immune to the dominant intramolecular dipole^ 
dipols mechanism and hence relax with a time constant much longer than 
the conventional 7\. The LPS states are generated by adiabatic transfer from 
non-equilibrium high-field states. The singlet states are allowed to decay in 
low field, and then transferred adiabatically back to high field where they can 
be detected as perturbations of conventional NMR signals. Low-field singlet 
states have been observed to decay with a time constant Tlfs — 104 ± 5 s, in 
a sample for which the measured 2\ is 16.8 db 0.3 s under identical conditions. 

The experimental procedure for investigating the lifetime of low-field sin- 
glet states is sketched in Fig-1. In its simplest form, the experiment applies 
to an ensemble of molecules, each containing a pair of spins-1/2 of the same 
asotopic type, but with different chemical shifts £i and i 2 , and a mutual J- 
coupling J. The experimental procedure uses a physical transport of the 
sample between two different magnetic field strengths, A,»gh and Sw, where 
■Shfch. » as shown in Fig.la. The transport processes from high to 

low field, and vice verses take times r^l^ and r^, Bp respectively, which are 
assumed to be less than Tj (both intervals are of the order of 10 seconds 



in the experiments reported here). The low-field storage interval tl.f, on 
the other hand, may be much longer than Ij, The transport between the 
two magnetic fields is synchronized with two radiofrequency pulse $equences 
denoted A4. and B in Fig. lb- These pulse sequences axe applied at the Lar- 
mor frequency of the spins in the high magnetic field, emd axe specified in 
Supporting Materiel 

The energy level picture at time point © in Fig-lc. shows the four Zee- 
man energy levels, -with energies in the order E aek < JS^ < < Eqq, 
assuming that 62 > Ji B the gyromagnetic ratio 7 and the ./-coupling axe 
both positive, and that the spin system is weakly-coupled in high field. In 
thermal equilibrium at ambient temperature, the lowermost level \ccot> has 
a small excess population (depicted by the filled balls), while the uppermost 
level 1/9/? > has a slightly depleted population (depicted by the white balls). 
These population differences are very small in the case of conventional ther- 
mal polarization, but may be of the order o£ one in spin hyperpolarization 
experiments (5-9). 

As discussed in Supporting Material, the pulse sequence A+ 1 which ifl given 
by two 90° pulses separated by a delay, exchanges the population of state 
|o**> with that of [&0> f and that of \J3a> with that of \fifl>. Hence, at 
time point 2, a non-equilbriuxn population distribution is generated, in which 
level ]arj9> has an excess papulation, while |>Sa> has a depleted population, 
as shown for time point ® in Fig.lc. 
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During the interval r4iap t the sample is transported to a region of low 
magnetic field J3icw, ivhich is of fete order of milliTesla. In low field, it is 
convenient to disciuss the behaviour of the nuclear spins using the four eigen- 
Stetes of the pure J-coapling HMftiltonian: 

Kj ~ 2*Jh ■ h [1] 

The four J-coupling eigenstates i»ay be classified as the three components of 
a triplet state, plus a singlet state. The triplet state components are denoted 

(rj> — \aa> 

[7h> = ^<|a0> + lfl*>) 
\T- 7 > - \Pfi> [2] 

and have the same J-conpliag eigenvalue of '+ Jar J. The singlet state 5s de- 
noted 

and has a J-couplIng eigenvalue of —^irJ. The triplet and singlet states are 
energy eigenstates of the nuclear spins only if the magnetic field is zero. The 
energy level rtiflgra™ i& this case is depicted next to time point ©in Fig-.lc 
Although we use the zero-field eigenstates as a basis for a representation of 



the spin density operator, the treatment below does not assume a field of 
exactly zero* 

If the transport from high to low field is fast compared to 21 , but slow 
compared to the /-coupling, the excess population of the high-field eigen- 
state \a&> is transferred fcriiabaticaUy into an excess population of the low- 
field eigenstate \Sq>. Similarly, the depleted population of the high-field 
dgenstate \0a> is transferred adiabatncally into a depleted population of 
the aero-field eigenstate \Tq >. Tfcese adiabatic transfer processes have been 
exploited before in the context of parahydrogen -enhanced NMR(7). The 
idealized distribution of the singlet and triplet state populations is shown at 
time point ® in Fig.lc This assumes pure adiabatic transfer and neglects 
T\ relaxation during the transport interval. 

The nuclear spina evolve in low magnetic field under a combination of 
coherent effects and incoherent relaxation processes. The coherent part of 
the low-field evolution has been studied before (10 f 11) but the analysis of 
incoherent molecular motion in very low field NMR has been restricted to the 
motional averaging of interaction tensors and to the investigation of lineshape 
perturbations (13-1$)> Here we use a standard second-order perturbation 
treatment of the modulated dipole-dipole interaction (16) to analyze the spin- 
lattice relaxation in low field. 
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la general, the spin density operator p evolves accor ding to the Liouville- 
von Neumann equation (1) 

w 

where the low-field Liouvillian is given by 

ii^ = r L p - i^lf [S] 

and ^Kjl^ is the sruperoperaior for the commutator with the low-field spin 
Hamiltonian: 

3C LF = -7^(1 + - 7Bw(i + + [6] 

The thermal polarisation terms £2 f) may he ignored in low field. If the 
relaxation is dominated by the i * tram olecular dipole-dipole coupling mecha- 
nism, the relaxation superoperafcor may be evaluated by second-PTder time- 
dependent perturbation theory and is given by (1) 

Ti* BS - f° K pr> (-T)3RCt» D {0) dr [7] 
J— 00 

where MpnOO is the commutation superoperator for the dipole-dipole cou- 
pling Hamiltopian, evaluated at time i, and the overtax represents an ensem- 
ble average. Eq,[4] is general for any field strength low enough to satisfy the 
extreme narrowing condition for the dipole-dipole relaxation. 



It Is convenient to use the matrix representation of the Liouvilli&n super- 
operator Zlp in a space of orthonormal spin, operators (18). The following 
set of six operators generate a suitable basis: 

|i) = \s Q xs Q [ - ivttf+ti*g--i?i*-irtit) 

|2) = \TiXTj] « /fjy 

|s) - |r Q ><r n [ = ^(/f/| + if/r + /, + / 2 -+A"4 + ) 

|4) « puxTLil = /f/| 

|s) = 4(|5 0 ><r 0 | + |r fl ><56|) - 

|e) = ^(|So><T Q |-|To><ffol) = -^{ite-SH) 

where J* are single-spin shift operators and J^} are single-spin polariza- 
tion operators (1,2). The ket |l) describes the population operators of the 
Singlet state- Kets j2^ , |s) and [4) describe the papulation operators of the 
three triplet states. Kets |s) and Je) represent the x and ^-components of the 
coherence between the singlet state and the central component of the triplet. 
All of these operators commute with the angular momentum operator along 
the field axis, as does the low-field Hamiltonian 3-C LF . The operator basis in 
Eq.[8j is closed under the application of the low-field Liouvxllian L LF . 

The matrix representation of the low-field relaxation superoperator r LF 
may be evaluated by standard relaxation theory assuming pure intramolec- 
ular dipole^dipole relaxation driven by rigid isotropic random rotation with 
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a correlation time r t . Tha low-field Liouvfllism evaluates to 
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18] 

where the dipole-dipole coupling between the spins is given by & » ^{jj^fA^-y^hr^ 
and r is the internuclear distance- The low-field chemical shift frequency dif- 
ference is defined 

Ln zero field (tjg" = 0), the first row and column of Eq.[8] contain only zsros. 
This indicates that the singlet sfcete population is dyrjjwnically isolated in 
low field, and is conserved under the evolution process. The sero-fiald singlet 
State may be used as a long-term repository for nuclear spin order, protected 
from intramolecular dipole-dipole relaxation. 

If the magnetic field is not exactly zero* Bq.[8] indicates that the singlet 
state is relaxed indirectly by dynamic couplings to the central triplet state 
and to singlet-triplet coherence. 



The idealised behaviour of the spin ensemble after em interval tl P in low 
field is depicted by the diagram at time point <g) in Fig-lc If the field is 
sufficiently low, and 7Lf Is long enough, the triplet pop-unions equilibrate, 
while the singlet population is conserved. 

The singlet population may be read out by adiabatic transport of the 
sample back into high field during the interval This transforms the 

low-field singlet population into a population of the high-field |o/3> state, 
as depicted for time point © in Fig.lc- The population of the |a/?> state 
is converted into observable NMR signals by the pulse sequence B. As de- 
scribed in Supporting Material, this sequence, which is given by three rf 
pulses separated by two unequal delays, converts the \c*0> state population 
into antiphase NMR signals centred at the chemical shift &. NMR signals 
deriving from T± relaxation during the transport interval, on the other hand, 
only generate NMR signals at the chemical shift 5 a . The interesting NMR. 
Signals deriving from low-field singlet storage are therefore cleanly separated 
from signals with a trivial origin. 

Experimental results for a 17.7 mM solution of 2,3-dibromothiophene in 
DMSO-d B at 20°C axe shown in Fig.2. The proton spin system is defined 
by <5i = 7.11 ppm, = 7,72 ppm and J = 5.7Hz. A conventional proton 
NMR spectrum, obtained at 400 MHz in the field B hlsh & 9-4 T, is shown in 
Fig.2a and displays the two doublets of a typical AX spin system. A con- 
ventional inversion-recovery series (not shown), shows that the spin-l&ttice 
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relaxation time constants for the inequivalent sites are very similar in high 
field = 16.5 ± 0.2 s; T t [tk) « 17.1 db 0.2 s). 

Fig.2b shows the Bpectrum obtained after the ©inglei>storage sequence 
in Fig.l, using a storage time r^p =r 100 a, and a storage field B im = 1.8 ± 
0.5 mT, The spectrum in Fig*2b shows diagnostic antiphase signals at around 
& = 7.72 ppm, and on antiphase dispersion structure around S z = T.llppnu- 

We can prove the phenomenon of low-field singtei storage by replacing the 
sequence by the sequence As discussed in Supporting Material* the 
sequence A^, which is derived from A+ by changing one of the pulse phases, 
exchanges the population of state |or£> -with that of instead of with 

[ac>. The population of |ckjS> is therefore depleted at time point© rather 
than enhanced. Fig. 2c shows the expected change in sign of the diagnostic 
antiphase signals around 62 — 7.?2 ppm. The signals around 5i » 7.11 ppm 
axe eJso perturbed, but in a more complicated way* since these signals are 
superpositions of $<-*veral contributions (see Supporting Material). Since the 
total interval r^j^p H-tlf4- r^Lsp = 120 s exceeds Tj by a factor of seven, the 
spectra in Fig.2b and c prove that the memory of the nuclear spin system 
has been extended significantly beyond the Ti limit. 

Fig.3 Bhows fche measured decay of the singlet population in low field, 
derived by monitoring the 7.72 ppm antiphase signals as a function of rjjp. 
In a field J3w = l.S ± 0.5 mT, the time constant for the singlet decay is 
estimated to be Tlps = ±53, The singlet decay time constant had 



no appreciable field-dependence for magnetic fields less than Bi ov , = 20 mT. 
This behaviour is expected from numerical calculations based on Eq. [8] (see 
Supporting Materia!). 

The singlet decay is observed to accelerate when the proton concentration 
in the solution is increased. This indicates that interrnolecular proton-proton 
relaxation is a strong contributor to the singlet relaxation, and that very 
long lifetimes may only be achievable in highly dilute solutions. Singlet 
relaxation is also expected to be caused by paramagnetic impurities such 
as dissolved ooqygen, interrnolecular interactions with solvent deuterons, and 
scalar relaxation of the second kind (16) via the Br nuclei. 

The long-term persistence of the aero-fleld singlet state cannot be at- 
tributed to the field-dependence of 21, which is very mild for the sample 
studied here. A more correct interpretation is that the low external field 
suppresses the chemical shift difference and renders the spins magnetically- 
equivalent, which isolates the singlet state from the triplet manifold, even 
with respect to dipole-dipole relaxation. La effect, the singlet nuclear spin 
state of 2,3-dibromothiopheae is an analogue of parahydxogen, which is also 
a nuclear spin singlet (6-8). 

The storage effect should be observable in systems of larger numbers of 
spins as long ee magnetic equivalence can be switched on and off by a change 
in the external magnetic field strength. We axe currently exploring fo\u> 
proton systems of the form X-CH2-CH2-Y which should display a similar 
property. 
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Singlet storage affects should be observable entirely in high magnetic 
field if a suitable rEMiiofreqiienCy pulse sequence is used to suppress chem- 
ical shift differences (19). Analogous !2 sero-field experiments in high field" 
have been performed in solid-sfeate Nb£R(20 r £l). Radiofrequency pulse se- 
quences might be designed which ©How the storage of nuclear spin order For 
times much longer than 7a, without removing the sample from the high-field 
magnet. 

Extension of the spin memory time should be useful whenever nuclear 
spin states are used as a long-term label for the position or chemical state 
of molecules. The long lifetime of lew-field singlet states is expected to 
facilitate nuclear hyperpolariz^tion experiments, especially those involving 
parahydrogen and may create new possibilities in the eacploitation of 

NMR for quantum computation (3,4)* 
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Figure Captions* 



Fig.l. (a) Sequence of magnetic fields u$ed.in the experiment. The sample 
is transported from a region of high, magnetic field -Su e u into a region of 
low magnetic field J3 1ow in a time r^ nBp . After remaining in the low field 
for a time tlf, ^ & sample is trsua$ported back into the higti-field region, 
a process which takes a time r£!Lp. 

(b) Radio-frequency pulse sequence performed at the Larmor frequency 
of the spins in high field. The pulse sequences and B are specified in 
Supplementary Materiel 

(c) Idealised sequence of spin state populations during the experiment. 
The four high-field energy levels belong to the states \/30>, \a>0> t 
)ora> and ]0ot>^ reading in clockwise direction starting -with the 
upper state. 
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Fig.2. (a) Conventional *H NMR spectrum of the solution of 
2,3-dibrorflothiophene (inset) in DMSO-d B - 

(b) Spectrum generated by the sequence in Fig.l, using the preparation 
sequence A+ and a storage time r uy = 100 s. 

(c) Spectrum generated by the sequence in Fig.l, using the preparation 
sequence A- and a storage time tlf » 100 a. 

Expanded sections of the spectra from 7,76 to 7.6Sppm are shown for 
(b) and (c). All spectra -were obtained from a single signal acquisition. 
The spectra in (b) and (o) axe expanded vertically by a factor of 2 
compared to (a), 

Fig.3. Decay of the diagnostic antiphase signals as a function of time **1f- 
The top and lower half indicate the integrated amplitudes of the two 
components of the antiphase multiplet centred at 7.72 ppiru The 
symbols represent measurements at different storage fields B\ w > The 
vertical scale represents the integrated amplitude of the spectral peaks 
as a percentage of their amplitude in an ordinary one-pulse NMR 
experiment. The solid lines are fits to decaying exponential functions 
with time constant Tjj?s = 104 s. The singlet decay is independent of 
the storage field in this regime. 
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Radiofrequency Pulse Sequences 
Sequence A+. 

The pulse sequence A+ is given by the two-pulse sequence 90<> — n — 90_oo, 
where the symbol 06 denotes a strong, non-selective pulse with flip anh 
gle fj and phase (bath angles are specified in degrees, and the phases 
taken into account the sign of the precession and the radio-frequeroy mixing 
scheme (1*2)). The delay Is set to the value Ti == |7r/w^ ffh | 4 where 

For the experimental results, the delay was set to the value tj 2,05 ms. 

Standard spin operator theory (3,4) shows tha-t in the case £% > 5i and 
7 > 0, tlie pulse sequence A+ has the approximate propagator 

U{A+) - «cjK-iJ(/i* + *.)} Hrf4*rJ*) [2] 

neglecting the effect of the ./-coupling during the interval Ti . This assumption 
is v alid in the case that the chemical shift frequency difference is much larger 
th£u> the J-coupling (weak coupling approximation). The transformations of 
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the individual state kets are 

U(A+)\accc> = -[(*&> 

U(A+)]0ct> = -i|c*a> 

U(A+)\00> - l£a> [31 

which leads to the exchanges of populations sketched in Fig.lc. The popula- 
tion of the \a&> state at time point ®is derived from the thermal equilibrium 
population of the \aa> state- 
Sequence A. 

The pulse sequence A_ is given by the two-pulse sequence 90iso ^ n — 90-bo- 
The sequence differs from A+. by a 180 p phase shift of the first pulse- The 
spin propagator for this sequence, under the weak-coupling approximation, 
is 

U(A^) = exp{+i|{/Lr + 1**)} exp{+ur/ lv } [4] 
which leads to the following ket transformations: 

U(AJ)\aa> = -\a£> 

U(A-)\fia> ~ i\0P> 

U(AJ)\tx0> hp i\ota> 

U(A-)\00> = \0a> . [5J 



^3 



The population of the |c*/? > state at time point (2>is derived from the thermal 
equilibrium population of the \&& > state. 

Since the \a{3> state is connected adiabatic&lly to the low-field singlet, 
the sequences A+ and A_ may be used to prepare low-field singlet states at 
time point (§) with populations deviating from the mean in opposite senses. 

Sequence B 

The pulse sequence B is given by 90 0 — — 1S0 90 — r a — 9O45, with pulse 
sequence delays given by r 2 - |ff/(2w£* h }| + 1 1/(47) | and 73 - |1/(4J)|. For 
the experimental results, the delays were set to the values t& = 44.64ms and 
73 = 43.61 ms. This sequence is designed to separate cleanly the interesting 
NMR. signals deriving from the low-field singlet state from the uninteresting 
' signals generated by relaxation of the spin system during the second transport 
interval. 

Rom standard spin-operator theory (3,4), the propagator U(B) for this 
pulse sequence transforms the population operator of state \a/3 > into a su- 
perposition of many operator terms, including some that represent antiphase 
single-quantum coherences: 

+ i «p{-w/4} Jfl? + J Hxp^iTr/4}/? Zf + . . . 

[6] 



The (— l)-quantum coherences represented by operators and gen- 
erate the diagnostic antiphase signals near the chemical shift 5 2 - 

Sequence B gives rise to no net signals at the chemical shift 5a when ap- 
plied to in-phase Zeenxan magnetisation, created by Ti relaxation during the 
second transport interval r^^. This may be seen from the transformations 

U{B)I U U{BY = -i«p{+i7r/4}j^/ 2 ft -iexp{^i7r/4}7 a -/? 

+| exp{+i*/4} I+Ii + \ *xpi-i*/4}Itli 

[7] 

and 

Partially-relaxed Zeeman magnetization of spin h only generates signals at 
the chemical shift while partially-relaxed Zceman magnetisation of spin T 2 
only generates multiple-quantum coherences, which do not generate a signal 
at ail In principle, croSS-correlated relaxation (5) during the transport in- 
terval could aim generate spurious signals at the shift Such signals would 
not change sign when exchanging the sequences A+ and A- y and appear to 
be negligible in the system studied here. 



Sample 

The samp le was 500 microlitres of a 17, 7 mM solution of 2,3-dibromothiophene 
in DMSO-d B , contained in a 5 mm high-resolution NMR tube equipped with 
a Young valve to facilitate degassing. The sample contained a water impurity 
in approximately 27 mM concentration. The sample was subjected to three 
freeze-pump-thaw degassing cycles, each lasting over 30 minutes, to remove 
dissolved oxygen. 

Experimental Procedure 

The NMR experiments were performed on a V&riaji Infinity^ 400 MHz NMR 
system at a field of B b igh ^ 9-4 T, using an actively-shielded 89 mm bore 
magnet. A standard 5 mm high-resohition NMR probe used. The 90° 
pulse duration was around 5 ^s. 

The sample was transported between the high and low fields using the 
standard pneumatic sample elevator. The air pressure of the elevator was 
increased during the lift operation, and reduced during the insert operation* 
in order to accelerate the transport times. With practice we could achieve 
fairly reproducible transport times of = 9 ± Is and rf®*Bp = 13 ± 2s, 

measured from when the sample is fully loaded into the coil to when the 
sample is outside the magnet bore. 



The field J9iow 'were estimated using a Hall-effect Gauss meter, with zero 
field calibrated to an accuracy of approximately 10 pT using a mu-metai 
shield far from the magnet. 

For a storage field 2?jaw = 20 roT, the sample was elevated to the top of 
the magnet bore and allowed to remain there for the time tlf- For the lower 
fields, the sample was carried physically to calibrated locations within the 
stray field of the magnet, or in the case 8ia V = 0, into the magnetically- 
shielded chamber. 

Numerical Simulations 

The Liouville-von Neumann equation may be integrated using the low field 
Liouv illian j&lp (Eq*[8] of the main manuscript) to predict the trajectory of 
the singlet population as a function of magnetic field and molecular motional 
parameters - 

For 2 a 3-dibromothiophene, the distance between the proton nuclei was 
estimated to be r = 253 pm. The spin-lattice relaxation in high field was as- 
sumed to be dominated by intramolecular dipole-dipole relaxation, modelled 
as the isotropic rotational diffusion of a rigid molecule. "With tins model, 
the observed T\ value of 20 b corresponds to a rotational correlation time of 
r fl «p 15 ps. 

Dynamic simulations of the singlet state population as a function of time 



and field are shown in Fig.l- At low values of the field Slow, the singlet state Fig-1 
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population ie predicted to be time-independent over a timescale of hundreds 
of seconds- As the storage field is increased beyond around 200 inT, a rapid 
low-amplitude oscillation of the singlet state population is superimposed on a 
slower decay. At higher fields the amplitude of the oscillations becomes very 
large and at long times the population equilibrates among the four available 
states, indicating the loss of all Spin order. 

The experimental results show a significant damping of singlet spin order 
on the timescale of 100 seconds for fields much less than 200 mT. This indi- 
cates the participation of relaxation mechanisms other than intramolecular 
dipole-dipole coupling. Intermolecular dipole-dipole relaxation involving wa, 
ter protons or protons on different DBT molewles is likely to he the major 
singlet damping mechanism. 
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Figure Captions. 

Fig.l- Simulated trajectories of the singlet state population for the case 
b\ = 7.11 pom, fa - 7.72 ppm, J = 5.7Hz, r = 253 pm and r c = 15ps. 
The values of Slow are shown in the figure- The population of the 
singlet state is tracked starting from an initial condition in which the 
singlet is 100% populated and all the other states have zero population. 
At high storage fields, the singlet state population is predicted to 
oscillate rapidly and decay towards a final value of 25%, representing an 
equal distribution of populations amongst the four states. The dark 
areas indicate very rapid oscillations of the singlet population. 
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